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1. Introduction 


The standard method of calculating the Bloch modes for a given lattice is to prescribe the 
constitutive relations for each point within the unit cell. The corresponding dispersion 
relations for the bulk material are then calculated either theoretically or numerically. 
One may interpret this dispersion relation in terms of an anisotropic permittivity 
and permeability ? Bu ik and //Bulk for a homogeneous medium. These, in general, are 
temporally and spatially dispersive. However, usually the permittivity and permeability 
e and jl of the materials used to construct the unit cell are considered to be only 
temporally dispersive, that is, they depend on the frequency lu [§} Up to now very little 
theoretical research has considered periodic media where the materials inside the unit 
cell itself are spatially dispersive, i.e. ?(cu, k) and k), where k is the wave number. 
A spatially dispersive medium responds not only to a signal at a particular point but 
to signals in the neighbourhood of that point. For this reason it is also known as a 
non-local medium. 

The use of a spatially dispersive inhomogeneous medium enables one to shape the 
propagation modes, extending the possibilities from the usual sinusoidal shape. This 
may have a number of advantages. For example in an accelerating drift tube, one could 
flatten the mode shape to enable higher acceleration for a given peak field. By contrast 
in signal transmission one may desire a higher peak for a given total energy. 

In order to construct an inhomogeneous unit cell it is necessary to consider a 
medium where the permittivity and permeability depend on both the wave vector k and 
the position x within the unit cell, i.e. 7(oj, k , x) and Jl(u ;, k , x). In order to make sense 
of the arguments of ?(cu, k , x) we no longer interpret ?(o;, k , x) as the Fourier transform 
of the permittivity response function e(t,x) but in terms of a differential equation [l]. 
This is given in equation (|3]) below. An alternative interpretation of ?(a;, k, x) is in 
terms of a susceptibility kernel [2,3]. 

This approach is in contrast with most theoretical articles which consider 
inhomogeneous spatially dispersive media. In 14-121 the inhomogeneity is restricted 
to looking at two homogeneous media which are connected and investigating the 
corresponding additional boundary conditions. 

For this article we will assume that there is a linear constitutive relationship between 
the polarization field P(t, x) = D(t, x) — toE(t, x) and electric field E(t, x). To simplify 
the analysis we assume: there is no magnetization so that H = ho'B and that all 
fields are functions of time t and one spatial coordinate x — x\, thus independent 
of x 2 ,x 3 . In frequency domain, k 2 = k 3 = 0 and we set k — k\. The electric and 
polarization fields are either linearly polarized transverse modes or longitudinal modes. 
For transverse modes E(t,x ) = E(t,x)e 2 , P(t,x ) = P(t,x)e 2 and B(t,x ) = B(t,x)e 3 
in the (ei, e 2l e^) frame. These assumptions automatically satisfies the two non-dynamic 
source free Maxwell’s equations V ■ B = 0 and V • D — 0 and the remaining Maxwell’s 


§ Note that in this paper we use ui to denote the temporal frequency rather than the angular frequency 
for which it is commonly used. 
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equations can be combined to give 


c 2 e 0 


d 2 E d 2 P d 2 E 
+ e 0 - 


dx 2 dt 2 ' " u dt 2 

Alternatively we assume that the electric and polarization helds longitudinal, and the 
magnetic held vanishes. I.e. E = E(t,x)e i, P = P(t,x)e i and B — 0 . Maxwell’s 
equations are automatically satisfied if 

cqE + P = 0 (2) 

The permittivity relationship between P(t,x ) and E(t,x ) analysed here is 


T , . -, . (3(u},x) 2 d 2 P - 

L(uj,x)P{uj,x) + — — (cu, x) = e 0 E(uj,x) 


(3) 


(27t) 2 dx 2 

where /3(co,x) is the finite propagation speed of the polarisation wave multiplied by 
a constant, giving it dimensions of length, and L(uj,x ) is a dimensionless constitutive 
quantity. The Fourier transform of P(t, x) with respect to t is given by 


P(co,x) = 


— 27 ricot 


P(t, x)dt 


(4) 


If L(oj,x) = L{pj) and /3(u, x) = (3(u>) are independent of x then we can Fourier 
transform pi) to give 


where 


P(u,k) = 


P(u,k) = 


e 0 E(u,k ) 
L(oj) — fd{uj) 2 k 2 


e - 2ni ^t+ k Pp(t, x ) dt d X 


(5) 


( 6 ) 


' —oo J —oo 


is the Fourier transform of P(t,x). Such a constitutive relation can be made 
using diagonally crossed wires 0- If L(oj) = ((a; — i\o) 2 + u}p)/ujp then © is the 
generalisation of a single-resonance Lorentz model of permittivity 13 to include spatial 
dispersion. For longitudinal modes, one can also use a wire medium 14,15 . Using the 
longitudinal component of? given by [l5j eqn. (4b)] [j 

-e 0 E(uj, k) 


P(u, k) = 


cu 2 /k p — k 2 /kl 


where, the grid spacing b in the same in both y and z and 
2 = 2ir/b 2 

P log b — log 27rr + 0.5275 


(7) 


( 8 ) 


For transverse modes, we consider a periodic structure where L satisfies 

L(uj, x + a) = L(oj, x) (9) 


|| Comparing notation: 27 r W ( here ) = fc(Song) and 2 lTk {he re) = ^(Song) We set a ( Son s) = &( Son g) = 
fe(here) 
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and /3(co,x) = /3(u>) is independent of x. We further assume that the amplitude of the 
inhomogeneity is dominated by the first mode, that is 

L(oj, x) = L 0 (lo) + 2A(cj) cos (27 t x/a) (10) 

In the case of longitudinal modes, it is the plasma frequency k-p given in (|8j) which is 
periodic in x, caused by a periodic variation in r. 

In both the transverse and longitudinal cases, the goal in this article is to find 
solutions to the source free Maxwell equations, i.e. the dispersion relations. However 
since the medium is inhomogeneous it is not possible to find single mode solutions of the 
form e 2m ^+ kx ) _ In section 1 we combine the polarisation differential equation ([ 3 J with 
Maxwell’s equation to get a single 4th order ODE. Since this has periodic coefficients, 
this is an example of Floquet’s equation. Hence we look for Floquet’s modes of the form 


P(t,x) = e 2niujt P(x) = e“ P <? 


0 2 tt i(q-\-K)x/a 


q =—00 


where n, 0 < n < 1 is the phase. This gives rise to a difference equation for P q 
A-(u))Pq- 1 + f q (u)P q + A (u)P q+ i = 0 
where in the transverse case 


f q (u) = L 0 (u>) ~ 


P{u) 2 {q + K) 


+ 


a 2 uj 2 


a 2 u 2 — c 2 (q + k) 2 


This is shown in the appendix. In the Longitudinal case, for a wire medium, 

r ( \ 2 ,2 (q+K) 2 

f q (uj) =u -k 0 -—— 

LL 

where kf } = k 2 | r=ro is a constant. 


( 11 ) 

( 12 ) 

(13) 

(14) 


We can trivially solve (12) by arbitrarily fixing Pq and P\ and then applying (12) 
iteratively to fold all subsequent P q . However, in general, this will lead to the coefficients 
P q diverging, i.e. \P q \ —> 00 as q —> ± 00 . The Fourier transform of this would therefore 
be a non-smooth wave which may not even be continuous. Therefore, for physical 
solutions, we demand that 

\P q \ — y 0 as q—t ±cxd (15) 


In section [3] we make the assumptions that the inhomogeneity A (u) is a small 
constant. This enables us to give a dispersion relation for uj and k, as well as the 


harmonics P q . In the unperturbed case A = 0, then (12) becomes f q (uj n )P q = 0, which 


we solve by setting uj n = D n and P q = S q where f n {D n ) = 0. We find that, analogous to 
non spatially dispersive media, there exist an infinite set of modes, (u n , P n (x)), where 
n G Z. Thus we write © 


as 


P n (t,x ) = e*™ nt F n (x) — e 


^2-iriUnt 


E 


q =—00 


P' ,v e 
r q V 


2 ni(q+K)x/a 


(16) 


where 


+ 0(A) 


( 17 ) 
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and 


/n(fin) = 0 


(18) 


We “normalise” (16) by setting 


P? = 1 


(19) 


The relationship between ui n and n which satisfy (18) corrected by the perturbation (17) 


is a dispersion relation. An example of the dispersion relations, both of the unperturbed 
system and the perturbed system, is given in figure [l] for the transverse case and figure 
[4] for the longitudinal case. 


We do not solve (12) exactly. Instead we find approximate solution to some order, 
that is we solve 


act- 1+/«K)cr+Acr+i=o(A r > 


( 20 ) 


for some order 0(A r ). In this article we solve (20) at least to order 0(A 3 ), which is 
where the new features appear. However in many cases we can do much better. For 
example, with the uncoupled modes given by ( J33| below, P r q l = 0(Al n_<? l) so that we 
can trivially solve (20) to order 0(Al n_,? l _1 ), although in fact we solve (20) to order 
0( A l"-«l+2). 

There are two classes of solutions which we call uncoupled modes and coupled 
modes. The uncoupled modes have the property that 


fq{VL n ) ± 0 for all q yl n 
and are given by 

oj n = Q n + 0(A) and 
Solving (fl2|),(fl3l).(Jl7]) and 


Pq = E + O(A) 


( 21 ) 


( 22 ) 


, gives rise to a dispersion relation between k and c u n . An 


example is given in figure [lj 

The existence of coupled modes is a new feature that arises with spatially dispersive 


media. These occur when (21) is not satisfied. That is there exist integers n and m and 


a frequency Q n G C such that n 7^ m and 
fni^ln) = fm{^n) = 0 

In this case the modes corresponding to P n and P m , (16), are coupled by the 


(23) 


perturbation. These are explored in section |4j As stated we wish to satisfy (20) to 
order 0(A 3 ). This can easily be achieved if n — m > 4, where we say that coupled 
modes have decoupled to order 0(A 3 ). However if n — m < 3 then new solutions 
appear. These require taking higher order expansions for co n and P q . Furthermore if 
n — m < 2 then the frequency is shifted. In the case when n = 0 one can use the 
symmetry in order to construct the odd and even modes. An example of the lowest two 
modes is given in figures [2] and [3] 

We see in figure [3] that the spatial profile of the lowest modes are significantly 
altered by the presence of the spatial dispersion. This suggests that with higher order 


harmonics cos(47nr/a) etc. in (10), one will be able to choose the electric field profile. 
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CU3 




Figure 1. The dispersion curve for the the transverse inodes for single pole resonance 
L(u) = (w — i\) 2 /ujp + 1 where a = 0.8, /3 = 0.4, A = 0.1 and uj p = 0.15. The lowest 
mode n = 0 (left) and several modes (right) The curves are the unperturbed frequency 
Re(flo), (solid black), Im(fl 0 ) (solid red), and the perturbed frequency with A = 0.75, 
Re(wo) (dashed black), Im(w 0 ) (dashed red). 



Figure 2. The odd P° (red) and even Pf (black) transverse modes for L 0 = 1, 
A = 0.75, a = 1, P = 1. In this case k = 0, wf = 0.399 and cof = 0.753*. 



Figure 3. The odd Pi(x) (black) and even —iP®{x) (red) transverse modes for Lq, 
A, a, /3, to® and ixf given in figure [ 2 ] 
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UJ 3 


C0 2 


CUi 


CO 0 


K, 


Figure 4. Dispersion curve for unperturbed (black) and perturbed (red) for 
longitudinal modes in a wire medium. Here a = 0.1, b = 0.02, tq = 0.001, and 
r i = r o /100. Only the n = 0 modes can we see the difference between the perturbed 
and unperturbed modes. 



Figure 5. Real (black) and Imaginary (red) longitudinal mode shape for n = 0 and 
k = 0.3. Here the values of a, b and ro are given in figure [4] but n = ro/5 to see the 
distortion. Here P° 2 = 0.027, P° x = 0.174, P 0 ° = 1, P[ J = 0.154, P 2 ° = 0.019. 


In section [5] we present an alternative numerical met 
frequencies and P(x). Finally in the concluding section [6 
inhomogeneous periodic structures. 


lod for finding the permitted 
we discuss the application of 


2. Floquet’s Modes 

For transverse modes, taking the Fourier transform of ([Tj) with respect to t gives 

( 2 v t)~ 2 E" = —lo 2 (c~ 2 E + p 0 P) ( 24 ) 

Here prime is the partial differentiation with respect to x. In most cases, in the following, 
we will not explicitly write the c o argument, in E, P, L and /3. 

Taking the Fourier transforms of (J 3 ]) and ( 24 ) with respect to x we get 

(u 2 c ~ 2 — k 2 )E{k) = —u 2 p 0 P{k) 
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Figure 6. Wire medium with periodically varying wires 


and 


k 2 p 2 P+(L*P)(k) =e 0 E(k ) 


Combining these into a single convolution equation gives 


00 


(L*P)(k)=[pV- : j r P ¥ j)P(k) 


(25) 


(26) 


or equivalently 


/ 3 2 d A P(x) + / + ^ d 2 P(x) + 2 dL(x) dP(x) 


(2ir) 2 dx 4 


c“ j dx 2 


dx dx 


/ d 2 L(x ) i 27tco 
V dx 2 


+ 


-L(x)P(x) + 


27rur 


P(x) = 0 


(27) 


in (t,x) space. We see that (27) is a version of Floquet’s equation. Hence P{x) is given 
by 


Ho = E p i 


0 27ri(q+K)x/a 


(28) 


q =—oo 


hence GD- In the appendix we show that from ( floj ), ( [26] ) and (|28]) the coefficients P q 
satisfy the difference equation (12) with f q given by (13). 


Observe that having higher harmonics, such as cos(47rx/a) in (10), will result in 


(12) being replaced by a higher order difference equation. Likewise letting fd depend 


periodically on x will also increase the order of the difference equation. 


2.1. Longitudinal modes in wire media with periodic variation 

The case for longitudinal modes is very similar. Using ([7]),([8]) with the radius of the 
wire varying periodically, see figure [6] 

^> = ^"K t .g + 2W(L/ a ) +a5275 ) < 29 > 

where k^ = kp\ r=ro . Substituting into (|8j) we see 

fcp = A)q + 2Acos(27t x/a) (30) 

We may consider simply setting r{x) = r 0 +ri cos( 27 T x/a), which will produce terms such 


as cos( 27 t x/a) 2 and will increase the degree of the difference equation (12). For small 
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enough rq then we can ignore the higher terms cos( 27 t x/a) 2 , then A « (rifcpfe 2 )/(27rro). 
Combining (25), (J2| and (28) we once again get the difference equation (12) where 
is given by (141). The dispersion relation for uj and k is plotted in figure El Observe 


the the perturbed modes are dramatically different near n = 0 and n = \. This is near 
the point where the modes couple. The detailed analysis of this behaviour is future 
work. An example distorted mode is given in figure [5j The result is a three dimensional 
periodic structure, with periods (a, b,b). This could also be analysed using numerical 
techniques, however this would be more time consuming. 


3. Approximate Analytic solutions 

As stated, we assume here that the inhomogeneity of the medium given by A in ( [Io| ) is 
small and give approximate analytic solutions to (12) which are valid to at least 0(A 3 ). 


Let be a solution to (18), / n (h2 n ) = 0. For each fl n we can find the corresponding 


Pg such that (12) is satisfied to some order given below. 


In the case of transverse modes, from (13) it is clear that f n ( oj) = F(Q,uj ) where 

P{u) 2 Q 2 


F(Q,u) = L 0 (uj) 


+ 


a 2 u 2 


(31) 


a* a 2 cu 2 — c 2 Q 2 

where Q — n + k and hence k = Q( mod 1) and n = |_QJ, the largest integer n < Q. 
Thus giving the shape of figure [1} The perturbation given in (32) below also has this 
property. In the longitudinal case F(Q,u) is replaced by F(Q,u >) = co 2 — — Q 2 /a 2 . 

In this section we consider the generic case, i.e. that there are no other integer m 
such that / m (Q n ) = 0. The special cases when there exists n ^ m G Z and D, n such 
that f n (I2 n ) = f m (D n ) = 0 are considered below in section |dj Let 

1 1 X + 0(A 4 ) 


UJn. — ^?7. + A. 


77 / 77 

J 77, n —1 


+ 


F' n Fn+l 


(32) 


and 


(—A) 


q-n\ 


n \q-n\ T 
k =1 ' r n+k 


g-nl+23 


P H = 

q 


(—A)l' ?_n l 

n \q-n\ 

k =1 ' r n—l 


+ 0( Al 


+ 0{ K\ q ~ n \ +2 ) 


q > n 

q = n 
q < n 


(33) 


where 


Fq = f q (n n ) , 


T’ = 
J q 


dfq 


du 


d 2 fq 


and T, = ^ 


(34) 


These (uj n , P™) solve the difference equation (12) to an order of A which depends on n 
and q as follows 

A P n q _ x + f q {u n )p; + AP g n +1 = 0(Al^l+ 2 ) for q ± n (35) 

and 

AP;- 1 + + A P„"+‘ = 0( A 4 ) (36) 
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The proof of (35) and (36) is given in the appendix. As in standard perturbation theory, 


improvement in the accuracy of (32) and (33) can be made by higher order corrections 


4. Coupled Modes 

In this section we consider the case when there are two distinct integers n and m and 


fi such that = / m (fi) = 0. By writing (13) as a quadratic for (q + ft) 2 we see that 

(q + k) 2 = C + (fl) 2 or (q + ft) 2 = C_(fi) 2 


where 




+ LoY _ 4fi 2 (L 0 + !) 




P 2 c 2 


(37) 


(38) 


One of two situation can occur. First, n and m correspond to the same root C±(H), so 
that 

(n + ft) 2 = (m + ft) 2 (39) 

This implies 2ft = — n — m and so either ft = 0 or ft = These will be considered in 
sections 14.51 and 14.61 below. 

In the second case set 


and 


(m + ft) 2 = C_(fl) 2 


(40) 


(n + ft) 2 = C + (fi) 2 
Then eliminating ft we have either 

(C+(fi)-C.(n))6Z or (C + (0) + C_(H)) G Z (41) 

Since, for most of cu, C + (co) ± C_(u;) are continuous function these situations will also 
occur in general. 


In the Longitudinal case, from (14) we see that only the case (39) can occur. 

Since two modes (u n , P™) and (cu m , P™) are coupled, we find that, in all cases, there 
are two new modes. 

It may be possible for specially chosen Lq{uj) and (3(oj) that there may exist three 
or four roots integer roots to ( Jl3| . However these will not be considered. Similarly 
one can construct L 0 (oj) and /3(u>) such that /^(H) = 0 or /^(f2) = 0. In which case 


formulae such as (45) will not be valid. These will also not be considered. 


In the special cases the dispersion curves become isolated points relating ui and ft. 
For example when ft = 0 and /(H) = 0 . 


Recall our goal is to solve (20) to at least order 0(A 3 ). We define the right hand 


side of (20) as Q ”, i.e. 

Qq — A-Pq-l + fq( U n)Pg + ^-Pq+1 
so that we require *2" — 0(A 3 ). 


( 42 ) 
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4-1. Decoupled coupled modes. 


In the case when n — rn > 4 then the coupled modes, decouple to order 0(A 3 ), then the 


two solutions have the same oo n given by (32) and (18). The corresponding modes are 
given by 

f {-A) q ~ n 


q—n+ 2 \ 


n nx* 


p n = 
q 


(-A) 


n—q 


n n-q -r 
k =1 ' r n-k 

0{ A n ~ m ) 

o(A n ~q - 2 ) 


+ 0(A 


+ 0(A n ~ q+1 ) 


q > n 

q = n 

m < q < n 

q — m 
q < m 


(43) 


This satisfies (42), with Q” = 0(A 4 ), Q" = 0(A n - m ), = 0(A n - m+1 ), Q" = 


0(A 9 n+2 ) for q > m, q ^ n and Q n q = 0(A n q 2 ) for q < m — 1. The proof is by direct 
substitution. Observe that in ([43]) there are terms that are given only to 0(A r ). These 


terms are always at least order 0(A 3 ) and furthermore it is clear that (20) is satisfied 
to at least order 0(A 3 ). 

The second solution, P™, is given by 

P™ — Pq’ where q’ = n + m- q (44) 


4-2. Coupled modes, n — m = 1 

In this case n — rn — 1 we have two solutions given by 

111 


— Al n + Auj' + A 2 


2P n _ 2 P, 7 l _ 1 


+ 


2 P 


n+1 


P' 


4 f T> V T> V 

J n—l^ n J n J n 


pH pn 

n _|_ ’ 


n—1 


n—1 


where 

J n = ±(KK-i )-' 12 

The corresponding modes are given by 


Pn- 3 = — A 
Pn -2 = A 
Pn -1 = “ 


Plu' 


P n—2pn —3 

px a 2 


+ 0( A s ), 


Pn +1 ~ 


Pn -2 

+ 

A 


P'n-\Pn-2 

A 


P" 


n —1 


P" 


4 PL 


+ 


p' 

•J n 


n —1 


n —1 


P' 

•J n 


n —1 


T" 

^ n —1 


4 PL 


^P'n + l< 


n —1 


P” 

J n 

4 PL 


4 P'n 2 P, 

p> 

*'n—l 


P' 

^ n 

2 P n -I 


+ 


P 7 


n—2 


+ 


P 7 


P 


n+1 


(P 


n+1 j 


+ 0(A 3 


2Pn+1 
Pn+2 = 


2Pn-2 

A 2 


Pn—2 
2\ 


Pn+2pn+1 


+ 0(A : 

+ 0(A 3 ) 


+ 0(A 3 ) 


(45) 


+ 0(A 3 


Pi = 1 , 


(46) 


Equations (45) and (46) satisfy (20) to order 0(A 3 ). This follow from direct substitution. 
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4-3. Coupled modes, n — m — 2 


Of all the cases the case when n — m = 2 is the most complicated. The two modes are 
coupled at order 0( A 2 ). Therefore it is necessary to take the second order perturbation 
to see the coupling. The solution for oj n is given by 


oj n — D n + A 2 


ol T 1 1 

+ 


77 / 77 77/77 

^ J n— 1 J n % ' 7i+l 

where a is the solution to the quadratic 


+ 0(A 3 


X*/ ^.2 i I Pn \Pn —2 , -T-! Pn—lPn x/ \ 77/ n 

+-2« + -x- + +-2-x-+ « “ + = 0 




n+1 


n —2 




n—3 


(47) 


(48) 


The mode is described by 
= A 2 


CD 


Pri—Spn —4 


D _ A " _A 2 x^+0(A 3 ), 


+ 0(A 3 ), -P n _ 3 — —A 


+i-3 


+—3 


Of -j - 1 


P n —2 = a + Acd 2 + 0(A 2 ), P n _i = —A+— - A 2 -?+ + 0(A 3 ), 


Pn ~ 1 ? -Pn+1 — 


P n—1 P n—1 

A +0(A 3 ), + +2 = A ~ +0(A 3 ) (49) 


-Tr* 1 \ / 7 ft-l -J-- -T-* 

J- n+1 n+2*'n+1 

We see that the A coefficients of P„_ 2 and the A 2 coefficients of P n _ 3 and P n ~i are 
related by the quantity co 2 . However it would require perturbing to order 0(A 4 ) to 
determine this quantity. 


4-4- Coupled modes, n — m = 3 

In the case n — m = 3 one solution is given by 

A 2 / 1. 1 

Wn — 1’n + VT 


77 / \ 77 77 

n V* 7 n—1 ^ ? 


n+1 


+ 0(A 3 


and the mode given by 

+- 4 = -A 2 


«3 

Pn —4 


, Pn—3 ~ Aq ?3 + 0(A 2 ) , 


+-2 = A 2 


P n — 1) ++1 — 


1 _ a 3 

Pn—lPn—2 Pn—2 

A 


+ 0(A 3 ), P n _i — 




+ 0(A ), ++2 — 


A 2 


(50) 


A 


P n -i 


+ 0(A 3 ), 


n+1 


Pn+lP n+2 


+ 0(A 3 ) (51) 


where a 3 cannot be determined order at order 0(A 3 ). 

The other solution is given by exchanging n and m, i.e. by (50) and (51) where 
m — n = 3. 


4-5. Coupled modes when k = 0 


When k = 0 then from (13) we see that f q (ui) = f- q {P) for all ui. Hence T q = T- q and 


T[ = P'_ . Thus since f n (£l n ) = 0 then f_ n (Q n ) = 0. If n ^ 1, —1 then we have the 


decoupled mode, with solution given by (43) and (44). 


























Spatially Dispersive Inhomogeneous Electromagnetic Media with Periodic Structure 13 


The remaining case is when n — 1. In this case some simplification arises and the 
solutions may be interpreted as odd and even modes. The frequencies are 

A 2 


L0 


-P — fii + 


+ 0( A 4 


wf-O i + ^(-l + !-')+0 (Al) 


(52) 

(53) 


and the corresponding P” 

Rl 


p 1,6 _ 

r q 


-2A/P 0 

R- q 


q > 0 
q — 0 
q < 0 


(54) 


and 


p i,o _ 

r q 


R\ 


-R\ 


-q 


q > 0 
q = 0 
q < 0 


(55) 


where R\ = 1 and 


pi (~A) 

9 ur=\^ 


q -1 


+ 0(A 9+1 ) for q > 1 


1 +k 


These satisfy Q} ±1 = 0(A 4 ), Ql = 0(A 3 ) and Q\ = 0(Al 9 l +1 ) for |g| > 2 and 


Q\ = 0( A 4 ). The proof following similar to (35) and (36) given in the appendix. 


This is consistent with (47) and (49) where from (48) a = ±1. 


4-6. Case n = \ 


When k= \ then from (13) we see that f q (u>) = for all uj. Hence IF q = 4F-\- q 

and T' q = F_ x _ q . Thus since ,/ n (O n ) = 0 then /_i_ n (H n ) = 0. 

Again if n > 3 then the modes decouple. However if n — 0 or n — 1 then the modes 
remain coupled and we can use the results of n — m — 1 or n — m = 3 respectively. 
Unlike the case when k — 0 no significant simplification takes place. However we note 


that (45) simplifies to 


cu r 


= A 

n -p/ 
^0 


A 2 


T" 

J n 


r.p 2 T'ppy 


+ 0( A 3 


(56) 


5. Numerical Approaches 


A numerical approximation scheme, which is valid if A ^ L 0 , is as follows: Choose an 
integer N > 2. Then assume that P m ~ 0 for \m\ > N thus truncating the infinite 
set of equation given by (12) to a set of 2 N + 1 linear equations for {P_ap • ■ ■, Pn }■ 
Write this in matrix language Mb = 0 where M is a (2 N + 1) x (2 N + 1) matrix with 
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1 

1 


■e 



Inhomogeneous periodic cross 
wire structure to create in¬ 
homogeneous periodic spatially 
dispersive medium in cavity. 
Protons 
Drift tubes 



Optimised electric field profile 


Figure 7. Drift tube accelerator with spatially dispersive media. The fields 
experienced by the protons are optimised. 


M k ,k = fk-N-i, M k)k _i = M fejfe+ 1 = A and b k = P k -N-i ■ Solve det(M) = 
values for oj. The corresponding null spaces give P m . 

For example, for IV = 2 the matrix equation is 


/ /(—2) A 0 0 0 \ 


( P ~ 2 ^ 


/o\ 

A /(-1) A 0 0 


P- 1 


0 

0 A /(0) A 0 


Po 

= 

0 

0 0 A /(1) A 


Pi 


0 

V 0 0 0 A /(2) ) 


V P 2 / 


w 


0 to obtain 


6. Conclusion 


In this article we describe a model for spatially dispersive media with a periodic 
structure. Two types of media are considered, one based on a single-resonance 
Lorentz model of permittivity, the other on a wire medium. Approximate analytic 
solutions to Maxwell’s equations are found for such media for small magnitudes in 
the inhomogeneities. We demonstrate the existence of coupled modes. Higher order 
correction should be possible by repeating the procedure. If one were to solve (20) to 
order 0(A 4 ) then one would expect coupled modes where n — m — 4 would remain 
coupled. 


By combining the solutions above with additional boundary conditions 4149,16 


one 


will be able to predict the behaviour or electromagnetic radiation as it passes though a 
slab of periodic spatially dispersive medium. 

It should be possible to construct a meta-material with such periodic spatially 
dispersive properties, especially the wire medium. These can be used to test that the 
packets given above are indeed supported. In addition, numerical simulations using, for 
example, FDTD codes in a periodic cell, can be made. 

As stated in the introduction, one of the key advantages for using inhomogeneous 
spatially dispersive periodic structures, is that one can choose the shape to the 
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propagating wave. For example, one may consider adding a spatially dispersive media to 
a drift tube in order to optimise the fields the particles experience. See figure [7J Another 
example, in data transmission, it to shape the wave to have a higher peak value. For 
example for longitudinal modes if E{x) = —P(x) = 3 -1 / 2 ( cos(27rx/a) + cos(67rx/a) + 
cos(107rx/a)), then the peak field increases by a/3 but the root mean square per length 
of E is unchanged when compared to the sinusoidal case. 
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Appendix A. 


Proof that the coefficients P q satisfy the difference equation (12): The Fourier trans¬ 
form of(28) 


P(k)= P,s(k-(q + 


(A.l) 


q =—oo 


Substituting the Fourier transform of (10) and (A.l) into the left hand side of (26) we 

have the convolution 
_ ^ /»00 ^ _ 

(. L*P)(k)= L{k — k')P{k')dk' 


J dk' ^ L 0 S(k - k') + A6^k-k'+^j+A6^k-k'-^j^ Y P q s(k' - 

e p < M* - tt) + H k - i± f ±i ) + A i k 


, Q + K 


q =—oo 

q + k - 1 


q =—oo 

by setting q' = q — 1 and q" = q + 1 


(. L*P)(k)= Y, P ^ 5 (k- q -^ 1 )+ Y P q "-iM(h- 


q" + K 


q =—oo 


X] P q , +1 A5(k- 


q' + k 


q ——oo 


relabelling q' and q" as q 


c L*p)(k)= J2 s ( k 


q + K 


( L^Pq + AP q —\ + AP q+ i) 


q =—oo 


Equating to the right hand side of (|26|) 

q + K' 

31 k — 

q =—oo 


°° I 

Y 5 { k ~ + Ai Vi + AiVi) 


hence 


L 0 - (3 2 k 2 + 


q =—oo 


LJ 




c o 2 — k 2 c 2 


q + K 


u 


Pq + A Pq-l + A Pn+1 — 0 


2 _ p c 2 1^1 1 ( ?~ i 1 9+ 1 


where k — (q + hi)/a. 


□ 
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Proof of (35) and (j36|. Expanding / n (u; n ) and using 


gives 


fnis^n) fnif^n) T {pJ n DjfjE n U n ) 


= A 2 


1 1 
+ 


Fn -1 


I , 


n+1 


2 

+ 0(A 4 J 


2 +' + 

^ n “ * * * 


Using this and (33) the left hand side of ( J36j ) becomes 

= a( ^ + 0(A 3 )) + A 2 (-?—+ -!—)+0(A 4 ) 

J n _i y- n+i' 


•F, 


n —1 


H 

' A 2 


-A 


•77 


n+1 


+ 0(A 3 


+ 


A 2 


T ' T 

n —1 J n+1 

= 0(A 4 ) 


+ 


A ' +^+o(A 4 ) 


T, 


n— 1 




n+1, 


hence (36). To prove (35) consider Erst the case when q > n, then using f q (tu n ) = 
T q + 0(A 2 ) we get 


A K-r + f,MP? + Ap t+1 = A A-i-.. x + 0(A«-” +1 ) 

\llfc=l -An+k 


T /<j+ n ) 


(-A) 


q-n 


n q-n -r 
k =1 •'n+k 


+ 0(A q ~ n+2 )) + A 


(—A) 9+1 ~ 

n q+l-n -r 
k =1 * r » 


+ 0(A 9_n+3 ) 


- ( — 1 " (a - /,(«„) A j + 0(A«-”+ 2 ) 

(a - (/„(!+) + 0(A 2 )) + + 0(A»-"+ 2 ) 


I q-l-n 
fc=l 

(~ A )*" 1 


n g-l-n -p 
k =1 ' r n+k 


— 1—n 


n g-l-n -r- 
fc=l y n+fc 


— 1—n 


(—A ) 9-1 

n q-l-n x- 
fc=l ' r n+k 

= 0(A q ~ n+2 ) 


(A - A) + 0(A q ~ n+2 ) 


The case for q < n follows likewise. 


□ 






















